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We derive simple necessary and sufficient conditions under which a quantum channel obtained 
from an arbitrary perturbation from the identity can be reversed on a given code to the lowest order 
in fidelity. We find the usual Knill-Laflamme conditions applied to a certain operator subspace 
which, for a generic perturbation, is generated by the Lindblad operators. For a weak interaction 
with an environment, the error space to be corrected is a subspace of that spanned by the interaction 
operators, selected by the environment's initial state. 



The ability to fully control physical systems so as to 
faithfully encode or transmit quantum information is of 
great interest both for fundamental research and applica- 
tions. The main obstacle is decoherence, a form of noise 
which rapidly destroys the quantum nature of a state. 
Quantum error correction (QEC) techniques consist in 
encoding quantum information into a physical system in 
such a way that it can be actively saved from such effects. 
This requires certain assumptions about the noise, one of 
which is generally that it is weak in some sense. In the 
most common framework, one assumes that arbitrary er- 
rors can happen independently on each qubit with a small 
probability, so that the likelihood of a combination of n 
errors happening at the same time goes as the nth power 
of that small probability. This intuition was generalized 
in Ref. [I] to any type of perturbative noise generated by 
a weak Hamiltonian interaction with an environment, or 
by a Markovian evolution. Nevertheless, our general un- 
derstanding of the correctability of a code still relies on 
the Knill-Laflamme (KL) conditions ^ which have been 
proven to be necessary and sufficient only when applied 
to the Kraus operators of an exactly known channel, and 
not in a perturbative framework. 

This problem was considered in Ref. ^ where suffi- 
cient conditions for perturbative QEC to lowest order in 
a small parameter were derived. Ref. [T] considered sim- 
pler sufficient conditions which amount to applying the 
KL conditions to a certain set of operators. Here we show 
that the KL conditions are necessary and sufficient when 
applied to a subset of these operators. 

Specifically, we assume that the noise is modeled by a 
quantum channel Af{p) = ^iP^i whose Kraus opera- 
tors Ni are convergent power series in a small parameter 
e, and such that N{p) = p for e = 0. Without loss of 
generality we write 

N,^5nl + eE, + ^F, + 0{€'). (1) 

Our result states that correction to lowest significant 
order in e is possible if and only if the operators 
\, El, . . . , Em, and hence their linear span, satisfy the 
KL conditions (Equ. [I2| ). Note that this excludes Eq. 
This operator must be anti-hcrmitian for the channel to 
be trace-preserving, and can be eliminated to order by 
the unitary e~'^^° . 



This result applies to any one-parameter family of 
channels expressed themselves (instead of their Kraus op- 
erators) as a power series in the parameter t. For the 
generic case, the first term to be corrected is of order t 
and has the Lindblad form. Our result says that the KL 
conditions must be satisfied for all Lindblad operators 
and the identity (Equ. |4]) . 

We also study the one-parameter family of channels 
resulting from a constant Hamiltonian interaction with 
an environment. In this case, the first nontrivial term 
to be corrected appears to order t^. It is sufficient for 
the span of the interaction operators to satisfy the KL 
condition as shown in Ref. , but knowledge of the initial 
state of the environment allows for the correction of a 
potentially smaller operator subspace. 

I. GENERAL PERTURBATION 

A quantum channel Af can always be written as 
-^(P) = J2i^iMtpMj where > and Tr (m/m^) = 
dij. Indeed a channel can be seen as a positive linear 
operator, the Choi matrix, whose eigenstates are linearly 
related to the operators Mi, with corresponding eigen- 
values Xi. Therefore, from linear perturbation theory we 
know that if Af is analytic in a parameter t, then so are 
the operators Mi and the eigenvalues A^. This implies 
that in general the Kraus operators Ni = ^/XiMi are 
given by power series in \/t. If M{p) = p for all p at 
t = Q (weak noise), the most general expansion satisfies 

N{p)^p + tC{p) + 0{e) (2) 

for 

C{p) = -t[H,p]+Y,L^pLl - -{L\u,p} (3) 

i=i 

and = H. This is the familiar term entering the 
Lindblad equation, although this equation holds also for 
non-markovian families of channel. It can be obtained 
by writing A^ and Mi as power series in t, with the re- 
quirements that the eigenvalue equation holds to zeroth 
order, and that the channel be normalized. In terms of 
the Kraus operators expansion (Equ. Il} this means that 
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e = "v/t, Eq = and Ei = Li for i = 1, . . . , m. Therefore, 
in this context our result states that necessary and suffi- 
cient conditions for correctability to order t are given by 
the KL conditions applied the all Lindblad operators Li 
(and the identity), namely, for all i and j, 

PL\LjP (X P and PL,P cx P. (4) 

II. WEAK INTERACTION 

Our result also applies to a model where the noise is 
caused by the system interacting with its environment via 
a Hamiltonian H which, for finite-dimensional systems, 
has the form 

n 

H ^\Y,'h®K, (5) 

i=l 

where Ji and Ki are Hermitian, and the operators Ji 
act on the system while the KiS act on the environment 
causing the noise. The parameter A > sets the overall 
strength of the interaction. 

If, in addition, we specify a certain initial state for 
the environment, we obtain an expression of the form 
given by Equ. [T]for the Kraus operators of the resulting 
channel. The terms of order e = t\ are 

n 

E,^iY,JMK,n (6) 

for i = 0, . . . , m, where — —1 and the states \i) form 
an orthonormal basis of the environment, including the 
initial state |0). If the initial state of the environment is 
mixed: p — '^iPi\4'i){4'i\i then one can use a purification 
|0) = Hi y/Pi\i^i) ® IV'i) and replace Kj by Kj (g) 1. 

In this context it was shown in Ref. [1] that for lowest- 
order correction it is sufficient to apply the KL condi- 
tions to the space Ji = Span(l, Ji, . . . , J„), which does 
not require the knowledge of the initial state of the en- 
vironment for its definition. Here we see that knowl- 
edge of the state |0) permits to identify the subspace 
Span(l, El, ... , Em) C J7i whose satisfaction of the KL 
conditions is both necessary and sufficient. In some cases 
where this subspace is smaller, this allows for better 
codes. An example will be given below. 



where \ip)p is any purification of p, ipp = |'0p)('0p|, and 
fip,a)^TT^Jp-2ap-2 ^um^\{^p\l^U\i;,)\ (8) 

is the fidelity between two states Note that some 
authors rather define the fidelity as the square of this 
quantity. Also, Fp{Af, id)^, where id is the identity chan- 
nel, is Schumacher's entanglement fidelity of J\f 5 . 

If the support of p is P, then Fp{J\f,A4) = 1 if and 
only if Af{a) = A4{a) for all <t supported inside P. 

We will also consider the worst-case fidelity, which does 
not depend on an input state p and is defined by 

i^r"(AA,7V/() := min Fp{U,M) (9) 
p=pP 

where the minimum is over all states p supported on P. 

The degree to which a code P is correctable can be 
measured by 

ap = maxFp(7^A/', id) (10) 

where P projects on the support of p, or alternatively by 
a'^'" = max i^^"" (7^A/', id) . (11) 

Clearly for e 0, A/" = id and both of these are equal 
to 1. We want to determine what their expansions in e 
are, and under what condition the first non-constant and 
non-zero term can be made to vanish. Note that although 
Fp{TZN , id) is linear in TZN , the optimal correction chan- 
nel TZ may also depend on e, making the whole expression 
non-linear in e. 

We find that the lowest non-constant order of both ap 
and ckp'" is e^, and that, in both cases, the terms of order 

vanish if and only if there exists \ij £ C such that 

PEjEjP = \,jP and PE,P = Aq^P (12) 

for alH 7^ and j ^ 0. 

This result is simpler than what one might have ex- 
pected, for instance considering Ref. [S] . In this work the 
authors attempted to find a modification of the KL con- 
ditions on the Kraus operators Ni. One can easily check 
that our conditions are equivalent to requesting the ex- 
istence of Xij € C such that, for all i,j ^ 0, 

PNlN,P^\o,eP + 0{e^) (13) 
PNjNjP = A,,e2p + o(£3)_ (14) 



III. MAIN RESULT 

We want to correct (reverse) the channel A/" on a code 
represented by the projector P, but we do not expect to 
be able to do it exactly. In order to quantify our success, 
we need a measure of how close two channels N and M 
are. We consider two such measures. For any state p one 
can define the entanglement fidelity 

Fp{N,M) :=/((AA®id)(^p),(X®id)(^p)) (7) 



IV. INDEPENDENT ERRORS 

Consider n independent copies of the channel A/" acting 
in parallel, then the overall channel on those n systems 
has Kraus operators 

n 

M^ = \^eY,E^^ +0{e^) (15) 

i=l 

Nhj^eEf +0{e^) (16) 
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where j ^ 0, and we write X^*^ for the operator X acting 
on the ith system. The other Kraus operators have order 
and can be neglected for lowest order correction. We 

see that all the operators Eg can be ignored since they 
only appear in the Kraus operator Mq which contains the 
identity component. 

For instance, we suppose that the noise on each qubit 
is due to an interaction of the form H = Jo(E)Kq + Ji(E)Ki 
with [Jo, Ji] 7^ 0. With no knowledge of the initial state 
of the environment, one is forced to use a quantum code, 
which may require 5 physical qubit to encode just one 
qubit. 

However, with knowledge of the environment's ini- 
tial state (say |0)), and our present result, we see 
that only Ei — « X)j matters. If fur- 

thermore E\ is normal (which happens here whenever 
Re(l|7^i|0)Im(l|ii:o|0) = Re(l|is:o|0)lm(l|ii:i|0)), then 
the linear span of 1 and E\ is equal to the span of 1 
and CTz defined in the eigenbasis of E\ . This implies that 
a simple repetition code can be used, necessitating only 
3 physical qubits to encode one logical qubit, or more 
generally any classical code correcting one error 



V. AMPLITUDE DAMPING CHANNEL 

We can readily apply our result to the amplitude 
damping channel on n qubits considered in Ref. 3 . The 
channel on each qubit has Kraus operators 



Nq = \-0[(^) and 7Vi=e|0)(l| 



(17) 



It describes the decay of a particle from an excited state 
|1) to a ground state |0) with probability proportional to 
. Given the analysis of the previous section we see that 
it is sufficient and necessary to correct the errors |0)(1| 
acting on any one qubit. However we have no simple way 
of reducing this to one Pauli error as |0)(1| is not normal. 



VI. PROOF OF THE MAIN RESULT 

It was shown in Ref. [7| that the quantity 



^P.- Tr , J2 NrP'N^Tr (aNjN,) (18) 



provides a good estimate for ap™ in the sense that 

min Ap,, < «™ < ^ min A,,, + ] (19) 

p=pP 4 p=pP 4 

for any state a — aP. 

This object Ap_o- measures how much information the 
"environment" obtains about the code. Indeed we can 
write (as in Ref. [7]) 



where A^(cr)Tr denotes the constant channel with output 
A/'(cr). The hats denote complementary channels; any 

channel M{p) = J^i-^iP^I can be written as result- 
ing from some unitary interaction with an extra system: 
the environment. The information that the environment 
then receives about the initial state is given by the com- 
plementary channel 



M{p) :=^Tr(M,pM;)|z)(j|. 



(21) 



Hence Ap^o- measures how close Af is from a constant 
channel, i.e. a channel carrying no information about its 
source. 



The fact that the expressions in Equ. [18] and Equ. 20 
are equal can be checked by using the definition of the 
fidelity as the right hand side of Equ. [8] and noting that it 
has the form max[/ |Tr {XU) \ for some operator X, which 
is reached when U is given by the polar decomposition 
of X, yielding maxu |Tr {XU)\ = Tr {VX^X). 

Equ. [19] implies that for any v > and any given a 
inside the code, Ap^o- = 1 + o{e'') for all p in the code if 
and only if a^'" = 1 + oie") . 

Similarly, it was proven in Ref. 8J that 



Ap,p < ap 



(22) 



Hence ap = 1 -f o{e'^) if and only if Ap^p = 1 + o{e'^) 

Our strategy is to compute the first nontrivial order of 
Ap.CT and find conditions under which this order vanishes 
in the case p = a. This solves the problem concerning 
Op. For ap, we will just note that the same conditions 
are clearly necessary for that same order to vanish, but 
also sufhcient, as a direct calculation shows they make 
Ap_CT independent from a to that order. 

We use degenerate perturbation theory to find the 
eigenvalues of the operator '^ij Nip^N^-Ti {aNlNj) up 
to order (see Appendix [A]). Expanding the square root 
of each eigenvalues in powers of e and summing them we 
obtain 

Ap,, = 1 + ^ Re [Tr {pE])Ty {ctEj)] 



^'E^Tr[(p + a)i?li?,] 



Ap,, = Fp(A/',AA(fT)Tr) 



(20) 



-t- Tr l^P^ E,p^ E] Tr [a{E[)^ E'^] + o{e^) 

(23) 

where P^ = 1 P projects on the kernel of p, and 

El:= E,-Tt{pE,)1. (24) 

We first consider the case p = cr as explained 
above. Let us define the non-trace-preserving (nor 
even trace-decreasing) completely positive maps £{p) = 
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J:^^oE^pEl and £'{p) = E.^o Using these 
definitions we can write 

Ap,p - l-e^Tr£'ip)+e^FpiV^,£'{p)Tr)+Oie^) (25) 

where 'P^(p) = P^pP^. The expression for the third 
term is obtained formaiiy in the same way that Equ. [20] 
is obtained from Equ. [18] 

Therefore the conditions for the first order term in the 
fidelities to vanish is 

Trf'(p) = Fp(^,f'(p)Tr). (26) 

Since the fideUty above is just the overlap between two 
states, it is bounded by the product of the norm of these 
states: 

F,{V^£,£'{p)Trf < Tr [P^£:(p)]Tr [£'{p)]. (27) 

Hence for Equ. [26] to be satisfied wc must have 

Tt£'{p) <TtP-^£{p). (28) 

Direct calculations shows that this means 

Tr {P£'{p)) = Tr {P£'{p)P) < 0. (29) 

But then P£'{p)P — because it is a positive operator. 
Indeed, it is a sum of manifestly positive operators of the 
form XX^i; 

P£'{p)P = Y,{PE[^p){PE[^-p)^ ^ 0. (30) 

i 

This in turns implies that each of the operators PE'^y/p 
vanishes: 

peIVp = (31) 

for all i 7^ 0, or simply, recalling that P projects on the 
support of p, 

PE,P cx P (32) 

for all i 7^ 0. 



Together with Equ. [26j this also implies that 

i^p(^,f '(p)Tr)2 = Tr [P^£{p)] Tr [£' [p)] (33) 

Hence the fidelity is maximal. Since also both the states 
that the fidelity compares have the same norm, they must 
be equal. This implies that both CP maps are actually 
equal, namely 

P^(a)-£'(p)Tr(a) (34) 
for all states a supported on P. This means that for all 

T^{P^E,aE^^)^T^[E[p{Er)^]Tv{a) (35) 
which is equivalent to 

Tr {PE]P^E,Pa) = Tr {E'^p{Er)^)Tr (Per) (36) 

being true for any state a, which in turn implies the 
operator equation 

PE]P^E,P ^ Tr [E[p{E'^)^]P (37) 

for alH 7^ and j ^ 0. Since PEiP cx P this also implies 

PEjEjP oc P. (38) 

Hence we have shown the necessity of the conditions ex- 
pressed in Equ. [12] The sufficiency is straightforward, as 
both conditions together can be easily checked to imply 

the equality of P-'-f and £'{p)- 
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Appendix A: Perturbation theory 

We suppose that 

N, = d,ol + eE, + e^F, + e''G, + .... 



(Al) 



For the channel Af{p) = J2i ^iP^i to be trace-preserving 
we need in particular that -Eg = —Eq- Let us write the 
operator 



5 = 5]iV,p2^jTr(c77V/7V,) 



(A2) 



as a power series m e: 

S = SQ + eSi + t^S2 + 6^53 + e^Si + O(e^). (A3) 
Note that 

5o=p' and 5i-[So,p2] (A4) 

It is straighforward but tedious to write the operators Si 
in terms of the operators Ei, Fi and the next order Gi. 
We will not write them explicitly here. 

Let \i be the ith eigenvalue of S. We have 

(A5) 



Let us write 



with the corresponding eigenstates 

|A.)=5]e"|n,z). 

n 

As a reference, we will use the basis 

We want the eigenvalue equation 
^|A.) - A,|A,) 



(A6) 



(A7) 



(A8) 



(A9) 



to be satisfied to all orders in e. For e = we obtain the 
equation 



So\i)^\f\i) 



(AlO) 



Hence \i) is an eigenbasis of 5o. We assume without loss 
of generality that \^^^ — for i > d and nonzero for 
i < d. Let 



(All) 



i<d 



be the projector on the range of 5*0. 
We also introduce the "propagator" 

D, = {Xfh - So)-' 



defined to send the kernel of A^'^^'l — Sq to zero. Also, for 
i > d we write 

D = D,^~S^\ (A13) 
For order n > and all z, we obtain the equations 

n 

ixfh-So)\n,i) = Y,{S„r-\th)\n-m,i) (A14) 

m— 1 

Multiplying by this yields 

n-l 



m—1 



(n) 

Hence, in order to find A- we need to know the states 
|m, i) , m < n. 

In order to compute AJ , consider the projector Pi 
on the degenerate eigenspace corresponding to eigenvalue 
A-°\ so that 

P^\^) = V). (A16) 
From Equ. |A14| we have 

n 

PtM = ^ A(^™-A^)l)|n-™,z). (A17) 

m— 1 

Note that in our case PiSiPi — 0, hence {i\Si = {i\SiPj^. 
It follows that 

Xf^ = {^\S,\^) + {^\S,~X'i'^\l,^) 

= {^\S,\^) + {^\iS,~X'^^)P,^\l,^) (A18) 
= {i\S2\t) + mSiD,Si\t) 

where we used that A^ ^ = 0. 

We will see that for i < d we do not need to go to 
higher order. Therefore we now focus on the cases i > d. 
For these terms, A-^^ = 0. Indeed, noting that Pi = P-^, 
we have 



P^K2P^ = 



(A19) 
(A20) 



(A12) 



where we defined 

K2 ■.= SiDSi + S2. 
For the next order, Equ |A14| yields 

P\2,i) ^ DSi\l,t) + DS2\i) (A21) 

and 

-So\3,i} = Si\2,i) + S2\l,i) + {Ss - Af))|j) (A22) 

from which 

X^^ = {t\SiP\2,t) + {i\S2\l,t} + {i\S3\t) 

= {i\ SiDSi + ^2 |1, i) + {i\ S3 + S1DS2 \i} 
= {i\{SiDS, + S2)P\l,i) + {i\S3 + S,DS2\t) 
= {i\SiDK2 + S2DS, + S3\i). 

(A23) 
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As it turns out, again 



P-^KsP-^ = 



(A24) 



where we defined 



for all i,j>d and i ^ j. This means that the states |z) 
ioT i > d must be eigenstates of the operators P-^K^P-^, 
where 



Ki K3DS1 + K2DS2 + S1DS3 + S4 



and Xf "' are the corresponding eigenvalues. 

Since Sq — p^, we have A^^-* — pf; where pi are the 
eigenvalues of p and \i) its eigenstates. 
Putting everything together, we obtain 



(A30) 



K3 S1DK2 + S2DS1 + 83^ K2DS1 + S1DS2 + 5*3 

(A25) 

and therefore — 0. Using 

P|3,z) ^ DSi\2,i) + DS2\l,i) + DS3\i) (A26) 

and 

- 5o|4,z) = 5i|3,i) + 52|2,i) + 53|l,i) + (^4 - \f^m 

(A27) 

we finally obtain 

Af' = li\SiP\i,i) + {i\S2\2,i) + {i\S3\l,i) + (^|54N) 

= (ili^a^^la, i) + (i| ^li?52 + ^3 |1, + (i| ^4 + ^1^^3 
= {i\K3P\l, i) + (i| 54 + S1DS3 + i^2i?52 |i) 

= {i\K3DSi + i^2£'5'2 + S1DS3 + 54 

(A28) 

This is the first nonzero term. We still need to know 



i>d 



i<d 



E 

i<d 



(2) 



1 2\ 

2 Pt 



i>d 



(A31) 



where a direct calculation yields 

i(2) 



more about the state From Equ. A27 



(j|ii:3DS'i + K2DS2 + S1DS3 + S4 \i) = (A29) 



E ^ = E Tr (p4)Tr - - ^ Tr ((cr + p)eIe, 

)i<d j>0 i>0 

(A32) 

and 

E = Tr /E P^S^P^i?]P^Tr(a(i?DtS^) 

i>d y ■ij>0 

(A33) 

where E'^ = Ei - Ti {(jE,). 



